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Abstract 

We consider the maximal displacement of one dimensional branch- 
ing Brownian motion with (macroscopically) time varying profiles. For 
monotone decreasing variances, we show that the correction from lin- 
ear displacement is not logarithmic but rather proportional to T 1 / 3 . 
We conjecture that this is the worse case correction possible. 

1 Introduction and statement of results 

The classical branching Brownian motion (BBM) model in K. can be 
described probabilistically as follows. At time t — 0, one particle 
exists and is located at the origin. This particle starts performing 
Brownian motion, up to an exponentially distributed random time. 
At that time, the particle instantaneously splits into two independent 
particles, and those start afresh performing Brownian motion until 
their (independent) exponential clock rings. 

We introduce some notation. Let Aft denote the collection of par- 
ticles alive at time t, set N(t) = \Aft\, and for any particle v € Aft, 
let x v (s),s £ [0,t] denote the (Brownian) trajectory performed by the 
particle and its ancestors. N(t) is a continuous time branching process, 
and it is straightforward to verify that N(t)e~ t is a Martingale, which 
converges almost surely to a positive, finite random variable rioo. In 
particular, t~ x log N(t) converges almost surely to 1. 
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We will be interested in the location of the maximal particle, i.e. 
in the random variable 

M t = max x v (t) . 

v£]\ft 

As is well known, the distribution function F(x, t) = P{M t > a;) satis- 
fies the Kolmogorov-Petrovskii-Piskunov equation (also attributed to 
Fisher) 

dF ld 2 F 

See |M75j for a probabilistic interpretation of the KPP equation. 
In a seminal work, Bramson |Br78j showed among other facts that 

M t = m(t) + Op(l), m(t) = y/2t-^~logt, (1) 

in the sense that for any e there is a K e so that 

P{\M t -m(t)\ >K t )<e. 

In particular, Med(M t ) = m(t) + 0(1), where Med(M t ) denotes the 
median of M t . (In subsequent work |Br83j . Bramson also discusses 
convergence to a shifted traveling wave, but this is not the focus of the 
current work.) Analogues of (JXJ) also hold in the setup of discrete time 
branching random walks, see [ABR09) : for a recent convergence result 
for BRWs, see [AiTT] . 

The leading term in ([T]), linear in time, is a relatively straight- 
forward consequence of large deviations computations and the first 
and second moment methods; in particular, the coefficient \[2 would 
be the same if instead of BBM, one would consider the maximum of 
e* independent Brownian motions run for time t. On the other hand, 
the logarithmic correction term in (TTJ) is more subtle, and reflects the 
correlation structure of the BBM: for the maxima of independent BMs, 
the 3/2 multiplying the logarithmic term is replaced by 1/2. For a 
pedestrian introduction to these issues, see the lecture notes |Z12) . 

Our goal in this paper is to address situations in which the dif- 
fusivity of the Brownian motion changes in time, in a macroscopic 
scale. This is motivated in part by our earlier work [FZllj . in which 
we showed that the corrrection factor 3/2 multiplying the logarithmic 
term can be replaced by different, and eventually much larger, factors. 
This naturally leads to the question, whether larger-than-logarithmic 
corrections are possible. Our goal here is to answer this question in 
the affirmative. 

We turn to the description of the time inhomogeneous BBM that 
we consider. Fix T (eventually, large). We consider the BBM model 
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where at time t, all particles move independently as Brownian motions 
with variance a^(t) — cr 2 (t/T), and branch independently at rate 1. 
Here, a is a smooth, strictly decreasing function on [0, 1] with range 
in a compact subset of (0, oo), whose derivative is bounded above by a 
strictly negative constant. Define Aft, {%v{t)}vejV t andM t , t £ [0,T], as 
in the case of time homogeneous BBM. This model has been considered 
before in IDS88] . Our main result is the following. 

Theorem 1.1. With notation as above, we have that 

Med(M T ) = v„T - g a {T) , (2) 
where v a is defined in ([5]). and 

< liminf 9 -^± < hmsup ^ < oo (3) 

We emphasize that it is already known a-priori [FalO that {Mt — 
Med(AfT)} is a tight sequence; in fact, the tails estimates in |FalO] are 
strong enough to allow one to replace, both in the statement above 
and in Theorem ITU the median Med (Mr) by EM T . 



2 Proofs 

Before bringing the proof of Theorem ll.il we collect some preliminary 
information concerning the path of individual particles. With W. and 
W. denoting standard Brownian motions, let 

X t = f a T (s)dW s , t£[0.T}. 
Jo 

Let r(t) — L a^(s)ds. Clearly, X. has the same law as W T ( t) ■ The 
following is a standard adaptation of Schildcr's theorem S66, DZ98 , 
using the scaling properties of Brownian motion. 

Theorem 2.1 (Schilder). Define Z t = ±X t/T ,t £ [0,1]. Then Z t sat- 
isfies a large deviation principle in Cq[0, 1] of speed T and rate function 

l( f ) = J JoiMl ds > /effi[o,i], 

I oo , else 

Here, -Hi [0,1] is the space of absolutely continuous function on 
[0, 1] that vanish at 0, whose (almost everywhere defined) derivative is 
square- integrable. 

We now wish to define a barrier for the particle systems that is 
unlikely to be crossed. This barrier will also serve as a natural candi- 
date for a change of measure. Recall that at time t, with overwhelming 
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probability there are at most e t+ °^' particles alive in the system. Thus, 
it becomes unlikely that any particle crosses a boundary of the form 
Tf(-/T) if, at any time, 



Mf) := f : p^rds>t. 
Jo 2cr2 ( s ) 

This motivates the following lemma. 

Lemma 2.2. Assume a is strictly decreasing. Then the solution of 
the variational problem 

*V:=sup{/(l):J t (/)<t,te[0,l]} (4) 
exists, and the unique minimizing path is the function 

f(t) =V2 [ a{s)ds. (5) 
Jo 

In particular, 



v a = V2 / a(s)ds . 
Jo 



(6) 



Proof of Lemma \2.Sl We are going to prove that no other functions 
can do better than /. That is, if some absolutely continuous function 
g satisfies g(0) — and the constraint Jt(g) < t for all < t < 1, then 
.9(1) < 7(1) = v a . In fact, denote <j>{t) = J t {g) < t for < t < 1, and 
then (j)'{t) — f^Ty a.e.. We can write g 2 (l) as 

,9 2 (1) = = VWQvtfidty ■ 

Using Holder's inequality, we have 

g 2 (l) < 2 (J <f)'{t)a{t)dt \ (J a{t)dt \ = V2v a (J <f> '(t)a{t)dt 

Using integration by parts, the above is equal to 

V2v a U{l)o-(l) - J <f>(t)a'(t)dt 

Since 4>(t) < t and a'(t) < for all < t < 1, the above is less than or 
equal to 

V2v a ^er(l) - J ta'(t)dt \ = V2v a J a(t)dt = v 2 a ., 
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where we apply integration by parts in the first equality. This com- 
pletes the proof. □ 
Proof of Theorem We begin with the upper bound in ([2]). The 
first step is to show that in fact, no particle will be found significantly 
above Tf(t/T). 

Lemma 2.3. There exists C large enough such that, with 

A = {3te [0,T},veM u x v (t) > Tf(t/T) + ClogT} , 
it holds that 

P{A) ^t^oo . (7) 

Proof of Lemma UTSl Recall the process X. in Co[0, T], whose 
law we denote by Pq. Consider the change of measure with Radon- 
Nykodim derivative 

dPi, ( f* f'(s/T) 1 f (f'(s/T)f 

The process X. under P is the same as the process X +T/(-/T) under 
P\. Note that for any t < T, 

f^ dX ,^ + ^f x ,^m ds . (9) 

Jo o-{s/T) o-{t/T) T J o- 2 {s/T) 

We then have, with r = inf{t <T:X t > ClogT}, on the event r < T, 



o ^(a/T) " a(t/T) T J o*(s/T) 

V2Ciogr 

and therefore, with r' = inf{< < T : X t > Tf(t/T) + ClogT}, we 
have, for k < T, 

(dP 

< Sp ^l Te[fe _ 1)fc) ex P ^ _^__ T JJ. 

Define 

9 = inf{i < T : there is v £ Af T so that x v (t) > Tf(t/T) + ClogT} , 
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and Zk to be the number of particles z G Mk such that x v (t) < 
Tf(t/T) + ClogT for all * < k - 1 and x„(t) > Tf{t/T) + ClogT 
for some fc — 1 < t < k. Then, 

T 

P{6 <T)<^P(0e [fc - 1, k)) < P(Z k > 1) , 
fc=i 

and, using a first moment computation, we obtain 

V2ClogT A 



a(0) 



P{Z k > 1) < < e fe P (r' G [fc — 1, fc)) < exp 
Therefore, 

P(0<T)<Texp ^ ^^ + 1 J- 

This completes the proof of Lemma 12.31 □ 
We need one more technical estimate. 



Lemma 2.4. With X. and C as in Lemma ] 2. 31 there exists a constant 
C G (0, 1) so that 

e T P (X t < Tf(t/T)+ClogT,t G [0,T},X T > T/^-C'T 1 / 3 ) ^»0. 

(10) 



Proof of Lemma \2~4\ Fix C G (0,1). We apply a change of measure 
similar to the one used in Lemma 12.41 whose notation we continue to 
use. We deduce the existence of positive constants c\ , ci (independent 
of T) such that 

Po(X t < Tf(t/T) + C log T, i G [0,T],X T > T/(l) - C'T 1 / 3 ) 
< e - T e c i(c"r 1 / 3 +io g T) 

■E Po ^exp X s ds^ l Xr >-c'Ti/3lx t <o : t<T^ , 

where here we used that —a' is bounded below by a positive constant 
and o~ is bounded above. By representing X. as a time-changed Brow- 
nian motion, the lemma will follows (for a small enough C) if we can 
show that for any constant C3 there exists a C4 = 04(03) > indepen- 
dent of C G (0, 1) such that 

D := E ^exp jf B s dsJ l Br >_ c , T i/3l Bt < ,i<T^ < e^ 7 "' 3 , 

(11) 
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where {B t }t>o is a Brownian motion started at — ClogT. Note how- 
ever that 



D < E (^exp Lf jf \ Bs \d^ 1| Bt |<tVs j < e —^ 3 , 

where here P. is a Brownian motion started at and the last inequality 
is a consequence of known estimates for Brownian motion, see e.g. 
[BS961 Formula 1.1.8.7, pg. 141]. □ 

Remark The estimate in (ITT]) can also be derived probabilisti- 
cally. Here is a sketch. It is clearly enough to estimate the expectation 
on the event F := {Leb({s : > B s > -C5T 1 / 3 }) > 1/2}. But P(F) 
decays exponentially in T 1 / 3 . Some more details are provided in |Z12j . 

We have completed all steps required for the proof of the upper 
bound in Theorem 11.11 Due to the strong tightness result in [Fa 10] 
and Lemma T2.31 it is enough to show that 

P({M T > T/(l) - C'T 1 / 3 } n A C ) -> . 

This follows from the first moment method and Lemma 12.41 

We turn to the proof of the lower bound. Call a particle v £ Nt 
good if 

Tf(t/T) - T 1 /' 3 < x v {t) < Tf(t/T), for all t < T. 

Set 

Mi ^ ^ 1^ is a good particle • 



Lemma 2.5. There exists a constant C > such that 
P(M > 1) > e" 



_ CT l/3 



Proof of Lemma \2.5\ Recall the process X. in Co[0,T], whose law 
we denoted by Po, and the measure Pi defined by ([5]). We then calcu- 
late the first moment 

EMI E ^ ~^-v is a good particle 

= e T P (Tf(t/T) - T 1/3 < X t < Tf{t/T), for all t < T) 



E Po 



ex P ( _ J a ( s /T)^ Xs ) ~*-{- Tl/3 < x t<0, for all t<T} 



Jo °(s/T) 

Repeating the computation in (0), we conclude that 

EM > cxp f-c 6 T 1/3 ) P f-T 1/3 < X t < 0, for all t < t) 
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Since under Pq , X. is a time changed Brownian motion, we have that 

P (-T 1 / 3 < X t < 0, for all t < t) > e^ 1 '* 
for some C7 > 0. Hence, 



Tnl/3 

£M > e~ CsT 

for some c 8 > 0. 

We next derive an upper bound for the second moment EM 2 . By 
definition, 

1 ^ / ^ i u , f' are good particle' 

When v =/= v' , we let be the branching time of the last common 
ancestor of w and v' . Then, the paths {x v (s)}o< s <t and {x v >(s) — 

X v ' {t vv i ) }f, t 

, < s <t are independent. Applying a change of measure sim- 
ilar to that used in the computation of EM, we can bound above 



EM 2 < EM 

+ f T e 2T - l E 
Jo 



■ exp 



exp ^ 
' V2 



T v% \ 

, -^jj^dX s - T I 1{-TV3<X1<0, for all 0<s<T} 



dX 2 — (T —t)\ 1{_T 1 /3<XJ<TV3, for all t<s<T} 



dt, 



where X 1 and X 2 arc two i.i.d. copies of X_ (under the law P ). The 
above is equal to 



EM + / E 



exp 



V2 



dX s lf-Ti/a^XKo, for all 0<s<T} 



lo <T(8/T) 

' ex P [~ j t cr(s/T) dX *) 1 {-T 1 ' 3 <X^<T 1 / 3 , for all t<s<T} 



dt. 



On the event {-T 1 / 3 < X] < 0, for all < s < T}, using integra- 
tion by parts, one has 



-dX] 



X\ 1 f T xy(a/T) J ^TV3 



ds < 



/ a( S /Tr' s CT (1) T7 cr 2 (a/T) " a(l) ' 
Similarly, on the event {-T 1 / 3 < X 2 < T 1 / 3 , for all t < s < T}, using 



integration by parts, one has 



it s " <r(l) + a(t/T) Tj t a 

3T l/3 



^ V( 8 /T) 
'is/T) 



ds 



<7(1) 
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Therefore, EM is bounded above by EM plus 



1{_T1/3<X1<0, for all Q<s<T} l{-T I /3<.Y2<T 1 /3, for all t<s<T} 







eft. 



The latter integral is less than or equal to 

T e 4Tl/3 /^dt<e c ° Tl/3 



for some eg > 0. Hence, one can apply the second moment method: 

P(M > 1) = P{M > 0) > > e~ CloTl/3 . 

This completes the proof of the lemma by letting C = ciq > 0. □ 
By a direct first moment computation (or using the already proved 
upper bound ([3])), the minimum position of particles at time AT 1 / 3 is, 
with high probability, greater than or equal to —AT 1 / 3 for a constant 
A' depending on A. Choosing A > C (with C as in Lemma l2.5j) . 
and using the independence of the motion of descendents of different 
particles in Af AT1/3 , the lower bound in ^ follows. □ 



3 Discussion 

Our choice of considering strictly decreasing diffusivity is not acciden- 
tal: the computations in |FZllj hint that this should correspond to a 
worse-case situation. In fact, we conjecture the following. 

Conjecture 3.1. Let a be a smooth function from [0,1] to a com- 
pact subset o/(0,oo). Then, ^ holds with v a determined by ^ and 
| 9ct (T)| = 0(TV3). 

In a slightly more technical direction, it would of course be of inter- 
est, in the setup of Theorem 1 1.1[ to show that 5 (T (T)/T 1 / 3 converges 
as T -> oo, and to evaluate the limit. Our methods are not refined 
enough to allow for that. 

Finally, we mention that results for the homogenization of the KPP 
equation are available, see e.g. (Nil, NRRZ12 and references therein. 
In the terminology we employ here, those results correspond to fast 
varying, or microscopic, time inhomogenuities. 
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